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THE RECOMBINATION EQUATION FOR INTERVAL PARTITIONS 


MICHAEL BAAKE AND ELHAM SHAMSARA 

Abstract. The general deterministic recombination equation in continuous time is analysed 
for various lattices, with special emphasis on the lattice of interval (or ordered) partitions. 
Based on the recently constructed general solution for the lattice of all partitions, the 
corresponding solution for interval partitions is derived and analysed in detail. We focus our 
attention on the recursive structure of the solution and its decay rates, and also discuss the 
solution in the degenerate cases, where it comprises products of monomials with exponen¬ 
tially decaying factors. This can be understood via the Markov generator of the underlying 
partitioning process that was recently identified. We use interval partitions to gain insight 
into the structure of the solution, while our general framework works for arbitrary lattices. 


1. Introduction 

Recombination is an important mechanism in population genetics. It is one of the funda¬ 
mental processes that underlies the time evolution of a population of individuals under sexual 
reproduction; see 131131E] for general background and naiaiin] for important contributions. 
Recombination is a stochastic process that, in continuous time and in the deterministic limit 
of large populations, leads to a nonlinear ordinary differential equation (ODE) in the Banach 
space of finite measures on a locally compact product space for the types that characterise 
the individuals; see mm and references therein. The flow of this ODE generally acts on 
an infinite-dimensional space, but has a number of nice properties, including preservation of 
positivity under the forward flow as well as norm preservation in the positive sector. 

Moreover, the recombination ODE admits a reduction to a finite-dimensional nonlinear 
system of ODEs via a suitable ansatz in form of a finite convex combination of probability 
measures that derive from the initial condition by the application of a finite set of (also 
nonlinear) operators, which are known as recombinators. The explicit solution of this ODE 
system was recently achieved in [5], building on previous work [3 111 13 [IS] , in a recursive 
way. The solution formula derived there was explicit enough to establish the exponentially 
fast convergence of the solution to a unique equilibrium that only depends on the initial 
condition, as expected from the previously known cases [7] as well as the general theory [3; 
see also the indroduction of [5] and references therein for more. 

The mathematical setting was based on the lattice P(S') of all partitions of a finite set, say 
S' = {1,2,... , n}, and was not restricted to the biologically motivated case of two parents, but 
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also allowed for an arbitrary recombination scheme. This results in an interesting nonlinear 
ODE system in its own right, and adds a nontrivial example of a solvable case. Nevertheless, 
the treatment in [5] did not consider other lattices and thus did not explore the full structure 
of the solution in lattice-theoretic terms. It is the aim of the present paper to continue this 
development, in particular by employing the lattice of ordered or interval partitions. This 
has the advantage that complete closed expressions can be derived for sizes up to four sites, 
and with some effort even for five sites. This way, the intricate and rather complex nature of 
the recursive approach becomes more transparent, and a better understanding of the solution 
is possible. Also, we are able to shed more light on the non-generic degenerate cases with 
singularities in the recursion, which were only briefly looked at in [S] and had been excluded 
from consideration in all previous attempts [laiiu]. Finally, motivated by the underlying 
partitioning process from [SJ Sec. 6] and its Kolmogorov backward and forward equations, we 
derive a linear ODE system for the solution of the recombination ODE, thus explaining why 
the nonlinear ODE was solvable in the first place and how the degeneracies emerge. 

The paper is organised as follows. In Section [21 we introduce the mathematical concepts 
from lattice theory that we need, with some focus on our later applications (with one result 
of general interest being added as an Appendix). This is followed by a review of some of the 
known results for the recombination ODE in Section [3l adapted to our more general setting. 
This section also includes some new results on the decay rates of the system. Then, we 
consider the ODE for the lattice I(S') of interval partitions in some detail in Section [H This 
lattice allows to explore the general solution without the combinatorial complexity of the full 
lattice P(S'), and leads to a number of previously unknown properties. Special cases (with 
up to five sites) are written out in explicit form (and with possible applications in mind) in 
Sections [5] and [6l before we sketch further extensions and directions in Section [3 including 
the derivation of a linear ODE system for the recombination process. 

2. Preliminaries 

Let S' be a finite set, for instance S = {1, 2,... , n}, and let P(S) denote the lattice of all 
partitions of S, with order relation so that A =4 B stands for A being a refinement of B. 
Below, we use the shorthands 1 = {S} and 0 = {{!}, {2}, ... , {n}} for the maximal and the 
minimal element of P(S). We use A AB for the coarsest refinement and AV B for the finest 
coarsening of A and B. 

Let V C P(S) denote a sublattice. It is called decomposable when there is a partition 
S = S 1 LJS 2 of S into two non-empty subsets such that V C P(Si) x P(S 2 ). Moreover, V is 
called reducible when a subset U C S exists such that V is isomorphic with a sublattice of 
P([/). In this case, one may find a partition {Ai, A 2 ,..., A^} G P(5'), with r < n, such that 
V can be identified with a sublattice of P({1, 2,... , r}). 

Clearly, we are primarily interested in lattices V that are both indecomposable and ir¬ 
reducible, which implies that 0,1 G V. The perhaps most important cases are P(5') itself 
and 1(5'), the sublattice of ordered or interval partitions. The latter are the partitions 
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A = { 741 ,^ 2 ) ■ ■ ■) Ar} such that all parts Ai are contiguous subsets of S, so 

^ = {{1,.. {h + 1,.. • ,^ 2 },.. .,{kr-i + 1,... ,n}}. 

Another interesting lattice consists of all non-crossing partitions. 

Apart from our lattice V for the set S, we will need the induced lattices for all 0 7 ^ [/ C S'. 
They are dehned by 

Y{U) := MIt; MgV}, 

where A\u denotes the restriction of A to U. It consists of all non-empty sets of the form 
Ai<r\U as its parts. With this notation, one has V = V(S). We write |.A| for the number of 
parts of a partition A, and supp(^) := IJ1=i Ai for the underlying set that is partitioned. 

In our setting, we also have a non-negative function 

P(S) ^M^o, 

whose values q{A) for A G P(S) have the meaning of recombination rates in the underlying 
continuous dynamical system we want to study. Of particular relevance is the sublattice 

V G P(S) I q{A) > 0 ), 

which is the smallest sublattice of P(S') that contains all partitions with strictly positive rates. 
To guarantee that V is both irreducible and indecomposable, we will usually assume that it 
contains all relevant partitions with two parts. 

It will be vital to our later analysis that quantities defined on V(S') have consistent coun¬ 
terparts on all induced sublattices for non-empty subsets U C S. For the recombination rates, 
this is done via 

(1) AA) := E 

bgv(s) 

This ‘top down’ formula is consistent with the standard marginalisation of probability vectors 
for the top system to probability vectors for the subsystem defined by U ; see 0 Sec. 5] 
for a detailed account. 

One quantity that requires a different approach is the decay rate of a partition, which 
needs a ‘bottom up’ recursion. The decay rates V’ = a-re dehned recursively via those of 
the subsystems for non-empty U C S. One has mu 

( 2 ) = 

A^l 

which implies = 0 whenever \U\ = 1, together with the recursion 

|T| 

(3) A < A ) := Y . A ‘ a \ A ,) 

i=l 

for any partition A = {Ai, A 2 ,..., A|_ 4 |} G Y{U). In particular, one always has V’^(Q) = 0- 
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Functions such as g ox ^ can be viewed as elements of the Mobius algebra Mjg(V) for V 
over the field M. This is the free real vector space with basis {e^ | ^ € V}, where is the 
function given by £_^{13) ■= g. Then, multiplication is defined by £4 • % — ^a^Ab'i compare 
[21 Sec. IV.4.A] for details. More important to us is the real incidence algebra [21 Ch. IV. 1] 

Ak(V) := {/: V X V —^ M | f{A,I3) = 0 whenever A ^ Bj. 

Here, addition and multiplication with a real number are defined as usual, while multiplication 
f * g is defined via the (generally non-commutative) convolution 

{f*g){A,B) = Y, f{A,C)giC,B), 

A^C^B 

where we follow |2] to mark the summation variable by placing a dot underneath it. Note that 
this gives (/ * g) {A, H) = 0 for A ^ B. Clearly, 5{A, B) = 6^ jg defines an element 6 € A^(V) 
that is the neutral element for multiplication, so 5*f = f*6 = f for all / G A]jj(V). 

Recall that an element / € is a unit (or invertible) when some g G A]^(V) exists such 

that f*g = g*f = S. This happens if and only if f{A, A) A ^ .A G V, and the inverse 

is then unique, hence written as f~^. The left inverse property gives f~^{A,A) = 1/f{A,A) 
together with the recursion 

= 77^ E r'(AC)/(C,H) 

^ A^C^B 

for A B. A similar recursion follows from the right inverse property. In the Appendix, 
we add a classic (but less well-known) non-recursive formula [16], which is sometimes of 
advantage. An important unit is the zeta function (, defined by C(A,B) = 1 fox A ^ B and 
by ((A,B) = 0 otherwise. Its inverse, fx = is the Mobius function of the lattice, which 
plays a central role in many (if not most) lattice-theoretic calculations; compare [21 Prop. 4.6]. 
It is given by /r(A, A.) = 1 for all A G V together with 

^l{AB) = - Y = - E 

A^C-iB A<ChB 

for A -< B. It can often also be given in closed form, which will be used below. 


3. ReCOMBINATON EQUATIONS 


The ingredients for the formulation of the recombination ODE are a finite set S together 
with an irreducible and indecomposable lattice V of partitions of 5, a general product space 
X = with all Xi locally compact, the Banach space A4(X) of finite Borel measures 

on X with total variation norm ||.||, and the set of nonlinear recombinators Rj^ for A G V, 
defined by 


Ml 


Rj^{u)) = 






.UJ 


i=l 


1 
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for oj € M{X). On the positive cone the recombinators satisfy Rj\^Rq = for 

A,B G V. Note that one can define a matching order relation among them by saying that 
Rj^ ^ Rq precisely when Rj^R^ = Consequently, one has R^ ^ Rq if and only if ^ ^ ;S. 
The recombination ODE now reads [3 El [5] 

(4) =Y1 ) 

AgV 

with non-negative numbers q{A) = ^>'^(^) as introduced above. The forward flow preserves 
the positive cone of Ad{X) and, within it, also the norm of the initial condition, Wq say, 
see [5] and references therein for details. Since this is also the setting for the dynamics of 
recombination, we restrict our attention to probability measures on X (denoted hy V{X)) as 
initial conditions from now on. 

The key to solve Eq. (j3|), with cjq G 7^(X), is the decomposition 

(5) ujt = a^iA) Rj^iujo), 

^gv 

which effectively separates the dependence on the initial condition from the time evolution. 
Here, the coefficients constitute a family of probability vectors on V with initial condition 
a^iA) = 6{A, 1). As a result of [5l Lemma 2 and Prop. 7], the induced ODE for the coefficients 
then reads 

\B\ 

(6) at(A) = -i/jiDatiA) + ^ q{B) Waf\A\B^), 

A^B^l i=l 

where is the marginalised probability vector on the sublattice Y{U). Note that this vector 
has a clear meaning in terms of a (stochastic) partitioning process, as explained in detail in 
[5l Sec. 6]. Each satisfies an ODE of the same type as in Eq. ([6]) for the subsystem defined 
hy 0 ^ U C S; compare [H Sec. 5] for details on the important marginalisation consistency 
of our equations. A step-by-step repetition of the results from in our more general 

lattice-theoretic setting now leads to the following fundamental result. 

Theorem 1. The Cauchy problem defined by the ODE ([1]) with initial condition ujq G A4{X) 
has a unique solution. The corresponding flow in forward time preserves the space V{X) of 
probability measures on X. For any cOq G V[X), the solution is of the form (l5|), where the 
coefficients constitute probability vectors on V that are the unique solution of the nonlinear 
ODE system ([6]) with initial condition Oq(^) = 5{A, 1). □ 

The main progress of [5] in comparison to previous work is the insight that, for generic 
values of the recombination rates g, an explicit recursive formula for the solution can be 
derived. In our present setting, it can be formulated as follows. Let be the rate functions 
defined above, and set 

«tV) :=Y.e^{A,B)e- 

B^A 


(7) 


■ip^{B)t 
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for A € V(C/) and all non-empty U C S, where 6^ is an element of the incidence algebra 
A]g(V(C/)). Now, let 0^(1,1) = 1 together with 

( 8 ) 

for any ^ ^ 1, which is just a reflection of the initial conditions we need. Inserting Eq. ([7]) 
into the ODE Q and observing the marginalisation consistency then leads to the recursion 


(9) 


9^{A,B) 


= E 


g^(C) 


|C| 


2 = 1 


for all .4 ^ ^ 1, provided that no denominator vanishes and that the exponential functions 

in Eq. ([7|) are linearly independent. The entire analysis from [5] for the lattice P(5) remains 
valid for our irreducible and indecomposable lattice V, and leads to the following result. 


Theorem 2. Let S be a finite set, and V = V(S') an irreducible, indecomposable lattice of 
partitions of S. Assume that the decay rates fi^{A) are recursively defined for 0 U C S 
according to Eqs. m and (I3D, with the induced recombination rates on subsystems being 
given by Eq. (HD- Assume further that the decay rates fi^{A) with 4. € V are distinct. 

Then, the exponential ansatz dZD solves the Cauchy problem of Eq. ([6]) with initial condition 
a\^{A) = S{A, 1) for all 0 ^ U C S if and only if the coefficients 9^ are recursively determined 
by Eq. ([9]) together with 0^(1,1) = 1 and Eq. ([8]). 

Sketch of proof. Once again, this is a step-by-step transfer of the arguments for the lattice 
P(5) from [5] to V C P(5). Irreducibility and indecomposability of V ensure that the latter 
contains 0 and 1, which also implies the corresponding statement for Y{U) and all U under 
consideration. 

In particular, given 0 U 0 S, the distinctness of all with B G V(C/) follows 

from the assumed distinctness of the fi^{A) with 4 € V by an application of O Lemma 5], 
which holds in the generality we need here. This allows to use the ansatz ([7]) together with 
a comparison of coefficients to conclude as claimed, where all appearing denominators are 
non-zero as a consequence of the distinctness of the decay rates on all levels. □ 

The condition on the distinctness of the decay rates fi^{A) with 4 G V covers the generic 
case (both topologically and measure-theoretically). When looking at the excluded cases, 
one realises that some degeneracies are ‘harmless’ in the sense that the solution formula ([7]) 
remains valid (which also means that the selected exponential functions still suffice), while 
there are also ‘bad’ degeneracies that lead to singularities in the recursion and render this 
approach incomplete (due to the necessity of further functions beyond simple exponentials; see 
[5l Sec. 9 and Appendix]). Below, we shall discuss this phenomenon in detail for V = 1(5'), and 
later derive an alternative explanation in terms of Jordan matrices for linear ODE systems. 
At this stage, we recall [SJ Cor. 6] and state it for our more general situation as follows. 
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Corollary 1. Let S and V be as in Theorem^ Then, the generic recursive solution extends to 
all recombination rates with ^ € V, such that holds for all 0 U C S 

and all B eY{U). □ 

In [5], we also defined a linear reference system for V = P(5). The definition extends to 
any lattice V that contains 1 as 

(10) x(-4) = ^ q{^) with 1] := {C G V I ^ ^ C ^ 1}, 

where it is understood that we sum over the partitions of V, subject to the condition specihed, 
and analogously for the subsystems defined by the non-empty subsets U C S. In particular, 
one always has x(l) = ipiX)- Then, a simple calculation leads to the following result. 

Fact 1. If 1 £Y, the decay rates x of Eq. (uni) satisfy the relation 

xiD-xiA) = e{B). 

A4B-<1 

Moreover, if Y = I(S') and if q{A) > 0 holds for all A with two parts, one has 

xil)-x{A) > 0 


for all interval partitions .4 7 ^ 1 . □ 

Whenever \U\ ^ 3, we had and one has = X^{\) for any 7 ^ 0, while no 

general relation was derived in [5]. To formulate one, we say that a partition B splits a set 
V YU when the restriction B\y has at least two parts. It is then clear what it means to say 
that B splits one or more parts of another partition A. Note that B can split a part of A 
without satisfying B ^ A. 

This notation allows us to rewrite Eq. (IlOp for V = Y{S) as 

x{A) = e{B) = ^ q{B) = ^ ^ q{B), 

^ splits at least n^l 13 splits n 

one part of A parts of A 

where all partition sums run over Y. The analogous relation also holds for all Y{U) with 
0 A U C 5. Observe that the last identity implies the relation 

x^d) = E Y(a> = E E = E 

i^.4ev(u) i^.4ev(;7) bgv bgy 

B\jj=A B splits U 
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Now, 

we can rewrite 

Ip as 







1^1 


1^1 


Ml 




HA) = 


E 

x-'AX) ' 

= E 

E 

q{B) 


i=l 


i=l 


i=l 

B splits 







at least A 

‘•i 


\A\ 

/ 







II 

M 

Ee(8) 

+ 

E 

e{B) - 

■f 

E e{B) + 


i=l 

^ B splits 

B splits Ai and 

B splits 3 parts of 



only Ai 

Aj for one j ^ 

i 

Af including Ai 

(11) 

= E 

n Y ^(^) = 

Y*{AB)e{B) 

1 



n^l 

B splits n 


BeY 





parts of A 

where the second to last step follows from a simple counting argument and is the 

number of parts of A that are split by B. If we now compare the two calculations and recall 
that ^ 0 , the following result is immediate. 

Lemma 1. Within a given lattice V, the decay rates tp and x satisfy the relation 

'ipiA) = x{^) + '^{n-1) q{B). 

n^2 B splits n 

parts of A 

In particular, one has ip{A) ^ xA) for all ^ G V. □ 

Since V’(l) = x(l) whenever 1 € V, Lemma [T] together with Fact [T] provides an alternative 
way to express the difference — ip (A) in this case as 

ipiD-ipiA) = Y - ^(n- 1 ) Y ^(^)- 

A^B^l ri^2 B splits n 

parts of A 

This formula clearly shows that, in general, ipd) — ip{A) for AAl can have either sign. 

Let us now look more closely into the lattice V = 1(5'), which is simpler than P(5), but 
still displays most of the relevant general phenomena. 

4. Interval partitions 

This section deals with the case that V = 1(5), where 5 is a fixed finite set, for instance 
5 = {1,2,... ,n}. 

Lemma 2. Let 5 be as above, and assume that g^{C) = 0 for all C G 1(5) with \C\ > 2. 
Then, for any 0 7 ^ C/ C 5, one also has q^{B) = 0 for all B G I(C/) with \B\ > 2. Moreover, 
ip^{A) = x^{A) holds for all 0 ^ S and all A G 1(17). 

Proof. Since the induced (or marginalised) rates for B G 1(17) are given by 

Q^i^) = Y 

^ei(S) 

A\^=B 


THE RECOMBINATION EQUATION FOR INTERVAL PARTITIONS 


9 


the first claim is a simple consequence of the fact that, liB ^ I(C/) with \B\ > 2, any A € I{S) 
with ^1^ = B must have at least three parts, too, due to the structure of the interval partitions. 

When one considers the result of Lemma [U for the lattice V = 1(17), one finds that each 
term in the difference 'ip^{A) — x^(7l) is of the form g^{B) for some B € 1(17) with \B\ > 2, 
hence vanishes by the first part of our proof. □ 

For the generic choice of the recombination rates on 1(5'), our general solution has the form 

( 12 ) afiA) = 

with 9^ determined recursively as explained above. Let denote the Mobius function of the 
lattice 1(5), which is given by g,^{A, B) = (—l)l-^l“l®l for ^ ^ .S and /i‘^(7l, 0) = 0 otherwise. 
The definition of for 1(17) with 0 ^ U C 5 is analogous. These functions satisfy the 
following product identity. 


Fact 2. Let A,B,C E 1(5) with A ^ B ^ C. Then, one has 


ia^{A,B) = llfi^^{A\c^,B\c,), 

i=l 


where r = |C|, so C = {Ci, C 2 ,..., C^}. 

Proof. Since 7l ^ C, we certainly have |7l| = |7l.|c. |, and analogously for B. Now, one 

simply calculates 


Y{h^^{A\c^,B\c^) - 


l-^lcj l^lcj ^ |-4lcJ l^lcj 


i=l 


i=l 


= (-i)l-^l-l^l = ^i^{A,B), 

which establishes the product formula. □ 

With this preparation, one finds the following simplification of Eq. el- 


proposition 1. Under the assumptions of Lemma ^ where g(C) = 0 for all C E 1(5) with 
|C| > 2, the solution formula (fT^ simplifies to 

B)pA 

with CLS above in Eq. COD. Moreover, this formula for of holds for all remaining choices 
of the recombination rates. 

Proof. In view of Lemma [2l it remains to prove that 9^ = pL^. This will be done inductively 
again, so we will show that 9^ = pP holds for all non-empty U 0 S. 

When |17| = 1, one has 1(17) = {1} and 0^(1,1) = 1 = Likewise, when |I7| = 2, 

the corresponding lattice is 1(17) = {0,1}. It is clear that 0^(0,0) = 0^(1,1) = 1, and one 
easily checks that 9^{Q, 1) = —1 together with 9^{1, 0) = 0, so 9^ = pF in this case, too. 
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Now, let US assume our claim to be true for all finite sets U with cardinality |17| ^ n, and 
consider a set S with IS"! = n + 1. We know that 6^{A,B) = 0 whenever ^ ^ ;B, so let 
A,B ^ 1(5) with A ^ B and assume hrst that B ~< 1 (we will deal with the remaining case 
later). Now, our recursion together with the induction assumption and Lemma [2] gives us 


9^{A,B) 



Q^jC) 


\c\ 


i=l 


E 




Q^jC) 

X'^(l) - X®(^) 


|C| 

{A\c^,B\c^) 


i=l 


E 


/(C) 

x^(l) - X^{B) 




X'^(l) - X^ilS) 




where the second line follows because C -< 1 implies that all parts Q satisfy \Ci\ ^ n, while 
the hrst step in the last line is a consequence of Fact [2j 

Now, using the assumption and Lemma [21 the last sum evaluates as 


E /(C) = E//) - E = x/l)-x/S), 

wherefore we get 9^{A,B) = n^{A,B) for all ^ ^ ^ 1. Finally, since 0/1,1) = 1 is 

automatic, we consider an arbitrary A 1 and hnd 


0V,i) = - E ^/Ac) = - E 

yiW-<l 

from our previous argument. This hnally also gives 9^{A, 1) = 1) from the Mobius 

inversion formula. 

The hnal claim is clear, because the solution now extends to the non-generic cases as 
well. □ 


This result explicitly shows the compatibility of the general solution with the special case 
from [HE], which was also mentioned in [5]. What is more, it also shows that the linearity 
situation is slightly more general in the sense that one can still allow some recombination 
rates beyond the two parent case to be positive. Let us look into this point in more detail. 


5. Systems with up to four sites 

When we deal with 1(5) for any non-empty finite set 5 with up to 3 elements, we know 
that 

= X^ together with 9^ = and r]^ = . 

This means that the general solution is given by 

(13) af(^) = E/(AS)e-^"('')* 

B)pA 
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for the initial condition Oq (^) = 5{A, 1), irrespective of the values of the rates g. For \S\ = 1, 
this gives = 1 (we drop the upper index here), while |5| = 2, meaning I{S) = {1,0}, results 
in x(l) = i?(Q) thus in the coefficient functions 

a^(l) = and 0 ^( 0 ) = 1 - e"^®*. 

For S = {1,2,3}, one has 1 = (123), 0 = (1|2|3) and decay rates 


A 

X 

(123) 

£.(1|23) + ^i(12|3) + ^(1|2|3) 

(1|23) 

^(12|3) + ^>(1|2|3) 

(12|3) 

^(1|23) + ^>(1|2|3) 

(1|2|3) 

0 


The coefficient functions now read 

0 ,( 1 ) = 

ai(l|23) = _ g-x(i)t^ 

ai(12|3) = - e-^®*, 

a^(g) = 1 _ e-x(i|23)i _ e-x(i2|3)t ^ g-x(i)i. 

Here, one can easily check that they indeed form a probability vector. Note that formula 
m also holds for P(5), then with the Mobius function // and decay rate x of IF’('S'); see 
[5] for more. Note that the lattice P(S') contains one additional partition, namely (13|2), in 
comparison to 1(5). Up to this point, it is obvious that the probability vector also satisfies 
a linear ODE system — a point of view we will look at more closely in Section [7] below. 

For S = {1,2, 3,4}, things start to get more complicated. The lattice 1(5), in obvious 
notation, is given by 



(1234) 




\ 

(1|234) 

(12|34) 

(123|4) 

1 X 


X 1 

(1|2|34) 

(1|23|4) 

(12|3|4) 

\ 

1 

X 


(1|2|3|4) 

where an arrow points towards the next refined partition. An application of Lemma [1] to 
V = 1(5) gives the rate function as 


(14) 


^P{A) = x(^) + (^^(1|23|4) + e(Q))(i(A (12|34)), 
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with X as in Eq. (jlOp . An application of the recursion formula now results in the following 
expression for 9, which we show in (upper triangular) matrix form for convenience, 



(1|2|3|4) 

(12|3|4) 

(1|23|4) 

(1|2|34) 

(123|4) 

(12|34) 

(1|234) 

(1234) 

(1|2|3|4) 

1 

-1 

-1 

-1 

1 

X 

1 

— X 

(12|3|4) 

0 

1 

0 

0 

-1 

— X 

0 

X 

(1|23|4) 

0 

0 

1 

0 

-1 

0 

-1 

1 

(1|2|34) 

0 

0 

0 

1 

0 

— X 

-1 

X 

(123|4) 

0 

0 

0 

0 

1 

0 

0 

-1 

(12|34) 

0 

0 

0 

0 

0 

X 

0 

— X 

(1|234) 

0 

0 

0 

0 

0 

0 

1 

-1 

(1234) 

0 

0 

0 

0 

0 

0 

0 

1 


Here, one has 


X 


g(12|34) 

- iA(12|34) 


with 


V'(l) 


V^(12|34) 


^l(12|34)-(£.(l|23|4) + ^l(l|2|3|4)). 


This clearly shows that x can have either sign in general. When the denominator vanishes, 
we are outside the validity of the simple solution formula; see [H Ch. 9] for details. We shall 
come back to this case below. 

Before we expand on the solution, let us stick to the case that x 7 ^ 0, where 9 is invertible. 
The inverse function 77 is still surprisingly simple, and given by r]{A,B) = C{A,B) unless 
A = B = (12|34), in which case one has 


r7((12|34),(12|34)) = ^ 


7/;(1)-iA( 12|34) ^ ^ g(Q) + g(l|23|4) 
^(12|34) ^i(12|34) 


So, one sees that 9 = fi (and thus '!] = () happens precisely when '0 = X) which is equivalent 
with the condition ^>(0) = ^>( 1 | 23 | 4 ) = 0. 


Corollary 2. In the case of the lattice I({1, 2 , 3,4}), we get 'ip{A) = x(-^) for all A precisely 
when ^)(Q) = ^)(1|23|4) = 0. Then, one also has 0 = /i together with rj = Q, and the linear 
solution formula applies to all choices of the remaining recombination rates. □ 


Now, since 4 is the smallest system size where 9 can become non-invertible (which happens 
when X = 0), we can go one step back and write down the ODE system for the functions 
ai{A) on I{S) explicitly as 

(16) aM) = -maM) + E E 

where we have again dropped the upper index. This ODE follows from [5l Prop. 7] by restric¬ 
tion to the lattice 1(5) and the insertion of the explicit solution for the possible subsystems, 
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which have at most 3 sites and are thus given as above. Alternatively, one can use Eq. ([ 6 |) 
for a direct derivation. With [5l Lemma 8 ], see also Fact [3] below, the solution of Eq. (flGll is 


at{A) 

with the function 
(17) 


ao(-4)e-’^(i)* + q{C) 


Eo{a,P;t) 


t e if a = /?, 

('e-/3t _ e”"*^, otherwise. 


Note that Eq is continuous in the parameters, and symmetric under the exchange a 13. 

Now, observing that V'(l) — y>{A) = x(l) — x{A) for all A 7 ^ (12|34) by Eq. (fTH) and using 
the first identity of Fact [TJ the solution simplihes to 


aM) = aoi^) 

+ ^{A, (12|34)) £>(12|34) Eo(iA( 1), V’(12|34); t). 


with £»(12|34) = (V^(l) — V^(12|34)) + ^)(1|23|4) + ^>(0). The initial condition relevant to us is 
aQ(^) = 5(^,1). Together with the upper Mobius summation, = 5(^,1), 

the last equation then further simplifies to 


(18) aM) = + /x(A(12|34))(^(g) + ^i(l|23|4))Eo(iA(l),i/^(12|34);f). 

B)pA 

Here, when £»(1|23|4) + ^)(g) 7 ^ 0, the last factor becomes (via I’Hospital’s rule, 

in agreement with the definition of Eq) when ^)(12|34) —>■ ^)(1|23|4) + ^>(0), which means 
'0(12|34) —>■ V'(l). This shows the appearance of a new function in the expansion when we hit 
one of the essential singularities of the recursive solution described earlier. Note that the full 
solution can still be extracted from the 0 -matrix, if one starts with a generic case and then 
rewrites the contributions with a prefactor ±x in terms of the function Eq. 

In general, when ^)(1|23|4) = £»(Q) = 0, we are back in the situation of Corollary [2l and 
thus in the realm of the linear scheme. Note that no restrictions occur for the values of 
Q on (1|2|34) or (12|3|4), so that the linear scheme actually slightly extends beyond 12 ( 5 ')- 
In any case, one has ai{A) ^ 0 for all A G 1(5) and all f ^ 0, due to the preservation of 
positivity under the forward flow, together with X)Aei(S) ®t(^) ~ ^ ^ 0 , as a result 

of fJ-{A, (12|34)) = 0 and YIa^b E) = <5(0, H) together with ?/)(Q) = 0, which is 

nothing but a direct verification of the norm preservation mentioned in the Introduction. 


Corollary 3. The solution of the recombination equation dH on Ai(X), for S = {1,2,3,4} 
with the lattice V = 1(5) and with a probability measure oJq G ^{X) as initial condition, is 
given by 

^Ai^o)^ 

A&(S) 
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with the convex coefficients (^) of Eq. (fT8|) . 


□ 


6. Recombination with interval partitions for five sites 

Let us turn to the case of five sites, S = {1,2,3,4,5}. Here, as a result of Lemma [TJ the 
decay rates are related by 

V'M) = x(^) + 

B splits two 
parts of A 


as a partition A G 1(5) can have at most two splittable parts. In particular, ^^{A) = xA) 
for any A with at most one non-singleton part. For the other partitions, we get 


(19) 


A 

(12|345) 

(123|45) 

(1|23|45) 

(12|3|45) 

(12|34|5) 


f){A) - xA) 


^(1|23|45) + ^(1|234|5) + ^(1|2|3|45) + ^(1|2|34|5) + ^(1|23|4|5) + ^(0) 
0(12|34|5) + ^(1|234|5) + e(12|3|4|5) + e(l|23|4|5) + e(l|2|34|5) + ^(0) 
^(12|34|5) + ^(1|2|34|5) + ^(12|3|4|5) + ^(0) 

^(1|234|5) + ^(1|2|34|5) + ^(1|23|4|5) + ^(0) 

^(1|23|45) + e(l|23|4|5) + e(l|2|3|45) + ^(0) 


In particular, we get ip = X once again for a few more cases than those covered in Lemma [2l 
because no conditions emerge for the rates f 3 (l| 2 | 345 ), ^ 3 ( 12 | 3 | 45 ) and ^(123|4|5). 

Corollary 4. In the case of the lattice I({1, 2,3,4,5}), we get ip{A) = xA) for all A precisely 
when the eight recombination rates vanish that appear in Eq. (USD. Then, one also has 6 = p 
together with rj = (p, and the linear solution formula applies to all choices of the remaining 
recombination rates. □ 

The recursive formula for 6 leads to the matrix shown in Table [TJ Here, we used the 
shorthands 


= C(12|34|5), X 2 = C(12|3|45), Xg = C(l|23|45), X 4 = C(123|45) and X 5 = ^(12|345), 

where C(^) := £'("^ 0 /(V'(llsupp(. 4 ')) “ V'(^O) with A! being obtained from A by removing 
all singleton parts, so (12|34|5)' = (12|34) and so on. In these expressions, the rates on 
subsystems are the induced rates according to Eq. ([1]). It is a somewhat surprising feature 
that the final expressions take a relatively simple, systematic form only after exploiting the 
recursive structure in this way. 
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(1|2|3|4|5) 

(12|3|4|5) 

(1|23|4|5) 

(1|2|34|5) 

(1|2|3|45) 

(123|4|5) 

(12|34|5) 

(12|3|45) 

(1|234|5) 

(1|23|45) 

(1|2|345) 

(1234|5) 

(123|45) 

(12|345) 

(1|2345) 

(12345) 

(1|2|3|4|5) 

1 

-1 

-1 

-1 

-1 

1 

Xi 

*2 

1 

X 3 

1 

-Xl 

— X 4 

-X 5 

-X 3 X 3 +X 4 -X 2 

(12|3|4|5) 

0 

1 

0 

0 

0 

-1 

-Xl 

-X 2 

0 

0 

0 

Xl 

X 4 

X 5 

0 

X 2 -Xi-X 3 

(1|23|4|5) 

0 

0 

1 

0 

0 

-1 

0 

0 

-1 

-X 3 

0 

1 

X 4 

0 

X 3 

—*4 

(1|2|34|5) 

0 

0 

0 

1 

0 

0 

-Xi 

0 

-1 

0 

-1 

Xl 

0 

X 5 

1 

-X 3 

(1|2|3|45) 

0 

0 

0 

0 

1 

0 

0 

-X 2 

0 

-*3 

-1 

0 

X 4 

X 3 

X 3 

X 2 —X 4 —X 5 

(123|4|5) 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

-1 

—X 4 

0 

0 

X 4 

(12|34|5) 

0 

0 

0 

0 

0 

0 

Xl 

0 

0 

0 

0 

-Xl 

0 

-xs 

0 

Xs 

(12|3|45) 

0 

0 

0 

0 

0 

0 

0 

*2 

0 

0 

0 

0 

—X 4 

-X 5 

0 

XS+X 4 -X 2 

(1|234|5) 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

-1 

0 

0 

-1 

1 

(1|23|45) 

0 

0 

0 

0 

0 

0 

0 

0 

0 

X 3 

0 

0 

—X 4 

0 

-*3 

X 4 

(1|2|345) 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

-xs 

-1 

Xs 

(1234|5) 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

-1 

(123|45) 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

X 4 

0 

0 

—X 4 

(12|345) 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

X 5 

0 

-Xs 

(1|2345) 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

-1 

(12345) 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 


Table 1. The 0-matrix for the interval partitions of five sites. 


If • ... • X 5 7 ^ 0, the function 9 is invertible, and its inverse r] is given by 



(1|2|3|4|5) 

(12|3|4|5) 

(1|23|4|5) 

(1|2|34|5) 

(1|2|3|45) 

(123|4|5) 

(12|34|5) 

(12|3|45) 

(1|234|5) 

(1|23|45) 

(1|2|345) 

(1234|5) 

(123|45) 

(12|345) 

(1|2345) 

(12345) 

(1|2|3|4|5) 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 1 

1 

1 

(12|3|4|5) 

0 

1 

0 

0 

0 

1 

1 

1 

0 

0 

0 

1 

1 1 

0 

1 

(1|23|4|5) 

0 

0 

1 

0 

0 

1 

0 

0 

1 

1 

0 

1 

1 0 

1 

1 

(1|2|34|5) 

0 

0 

0 

1 

0 

0 

1 

0 

1 

0 

1 

1 

0 1 

1 

1 

(1|2|3|45) 

0 

0 

0 

0 

1 

0 

0 

1 

0 

1 

1 

0 

1 1 

1 

1 

(123|4|5) 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

1 

1 0 

0 

1 

(12|34|5) 

0 

0 

0 

0 

0 

0 

1 

Xl 

0 

0 

0 

0 

1 

0 J- 

Xi 

0 

1 

(12|3|45) 

0 

0 

0 

0 

0 

0 

0 

1 

^2 

0 

0 

0 

0 

1 1 

X 2 ^2 

0 

1 

(1|234|5) 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

1 

0 0 

1 

1 

(1|23|45) 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

^3 

0 

0 

i 0 

^3 

1 

1 

(1|2|345) 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 1 

1 

1 

(1234|5) 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 0 

0 

1 

(123|45) 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

i 0 

^4 

0 

1 

(12|345) 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 — 

^5 

0 

1 

(1|2345) 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 0 

1 

1 

(12345) 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 0 

0 

1 
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While the ^-coefficients suffice to write down the solution in the generic case, where we have 
'^^(1) 7^ foi' s-ll all non-empty U C S, the degenerate cases require some 

care. First, for A € {1, (1|2345), (1234|5), (1|234|5)}, the solution is always given by Eq. (fT3]l . 
These cases correspond to the four rows in the 0-matrix of Table [T] with constant entries, 
which agree with the values of the Mobius function. For A = (12|345) and A = (123|45), the 
general solution reads 

+ {il^iA) - x{-^))EQ{ip{l),'ip{A);t) 

with the function Eq from Eq. ([HD. Depending on the values of the recombination rates, one 
can thus pick up a term of the form as in the case of four sites. Next, one finds 

0^(1121345) = _ g-b(i|2345)t _ ^(i2|345) Eo(V’(l), V'(12|345); t), 

at(123|4|5) = - e-^(i234|5)t _ ^(i23|45) Eo(V’(l), iA(123|45); t), 

where the last term can be split as in the previous case. Similarly, one gets 
at(12|3|45) = ^)'(12|45)Eo(V’(l),V'(12|3|45);t) 

- g{l23\A5) Eo(iA(l), i/^(123|45); t) - ^(12|345) Eo(iA(l), V^(12|345); f) 

with £»'(12|45) = ij(12|3|45)-|-^(123|45)-|-^)(12|345). Here and below, the ' indicates an induced 
quantity on a subsystem (here with four sites). 

All remaining partitions are true refinements of (1|2345) or (1234|5), wherefore we may pick 
up a term of the form te“"* already in the ODE, depending on the choice of parameters. To 
deal with this complication, let us first state a simple variant of the results in [5l Appendix]. 

Fact 3. Let p as well as cri,...,a^ and be arbitrary non-negative numbers, and 

ip a continuous funetion on M^g sueh that eP^ip{t) is integrable. Then, the Cauchy problem 
defined by the ODE 

m n 

g = -P 9 + pit) s' 

i=l j=l 

together with the initial condition g{Q) = Qq and £i,£j G K has the unique solution 

/ f m n 

eP^p{s)ds -h - e"^*) -F ^ e'Eo(p, cr/; t), 

i=i j=i 

which holds for all t 0. □ 

Let us now consider the case A = (12|34|5) in more detail. According to Eq. Q, the ODE 
we have to solve here, after some calculations, reads 

= - HDaM + Y1 

+ cEo(V>(1234|5),iA(^);t) + ^ p{A,B){fi{B) - 

A^B^l 
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where c = ^)(1234|5) (^'(1|2|3|4) + ^>'(1|23|4)). This is of the type covered by Fact O with 
(fit) = c-Eq(V^(1234|5), t)- Note that the two rates here generically differ from p = V’(l)- 

Depending on the choice of the recombination rates, the solution can be a sum of exponentials 
(with coefficients given by 0), which is the generic case. However, it can also contain a 
contribution of the form (when x{^) 7^ V’('S) = '0(1) occurs) or even of the form 

which happens for the double degeneracy that = '0(1234|5) together with 

= ■0(1), as a consequence of O Lemma 9]. We leave it to the reader to analyse the 
details of the various cases that are possible here. Clearly, the treatment of the partition 
(1|23|45) is completely analogous, so that we have covered all partitions with up to three 
parts. 

The remaining five cases can be analysed in the same way. Doing so, one realises that the 
general form of the ODE for a^{A) with arbitrary A € I(S') reads 

aM)=-matiA)+ f^iA,B)[{^pa)-^P{B)) + {^PiB)-x{B))]e-^^‘^^^ 

+p(A(12|34|5)) cEo(0(1234|5),V’(12|34|5);t) 

+ /r(A(l|23|45))cEo(0(l|2345),V’(l|23|45);t), 

with c = £.(1234|5)(£.'(1|2|3|4) + £»'(1|23|4)) and c = ^(1|2345)(£»'(2|3|4|5) + £»'(2|34|5)). The 
solution is fully covered by Fact [3l and the necessary case distinctions should be clear. In 
particular, the solutions will contain simple exponentials, contributions of type Eq (whenever 
'i/j{B) 7 ^ xiB) in the sum, which leads to terms of the form te““*), but possibly also integrals 
as in the previous two cases (with the same consequences, namely possible terms of the form 
t 2 ^-at double degeneracies). Let us summarise this as follows. 

Corollary 5. The solution of the general recombination equation dl on Ai(X), for the set 
S = {1,2,3,4,5} with the lattice V = I(S') and with a probability measure ujq € V{X) as 
initial condition, is given by 

adA) Rj^iuJo), 

A&(S) 

where the convex coefficients a^{A) are the unique solution of the Cauchy problem defined by 
Eq. ([2011 together with the initial condition aQ(^) = <5(^1,1), as given in Fact El □ 

Beyond five sites, the complexity of the solution in the presence of degeneracies increases. 
In particular, for n sites, one can in principle also obtain terms of the form t^e”"* for any 
0 ^ m ^ n — 3, though this requires an m-fold degeneracy in analogy to the appearance of 
double degeneracies in our above examples. Since the generic case without any of the ‘bad’ 
degeneracies most likely covers the practically relevant cases, we do not want to expand on 
this issue any further. We shall come back to it from a different point of view shortly. 
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7. Further directions 

Let us return to the general case of n sites, with an irreducible and indecomposable parti¬ 
tion lattice V. At this stage, we know that (‘^) solves the recombination 

ODE (|3]) when each a-i{A) solves Eq. ([ 6 ]), both with the properly corresponding initial con¬ 
ditions. In the generic case with aQ(^) = 5{A, 1), this means a^{A) = 

where one also has C) (C)- This structure thus provides a mode de¬ 

composition (also known as decomposition into linkage disequilibria) = X^^tey (.A) with 

^,(A) = 

C^A 

= Y,v{AB)a,{B)Y,0{C,A)Rc{u;o), 

B)pA 

where the (signed) measure {A) satisfies the linear ODE 

i>t{A) = -il;{A)u^{A) 

together with the initial condition Ug(^) = X^c ^_4 ^(C,-A) i 2 g.(t(jQ). 

Since our general assumptions on the recombination rates (for V being indecomposable and 
irreducible) imply that 'ijj{A) > 0 for all ^ 7 ^ 0 , we get the norm convergence 

OJt 5" Rq{u}q), 

with exponential decay of all modes Vq{A) for ^ / 0. For any cjq E V^X), there is thus 
the unique equilibrium, R^^^ujq), and the convergence towards it is exponentially fast. In the 
degenerate cases, as explained above, this situation may have to be modified to include mode 
decay of the form for some m, where m ^ n — 3 with n the number of sites, but this 

does not change the equilibrium. As the general structure should be clear, we leave further 
details to the reader. 

Let us return to the generic case that 7 ^ ip^{A) holds for all A E V(C/) and all 

non-empty 1/ C S'. Here, we are in the regime of the recursive formula ([9]) for 9^. Though 
it seems difficult to write down a closed formula for 9^, we can say more on the diagonal 
elements. Define 'd^{A) = 9^{A,A) and set 'd^[0) = 1 . This way, d is an element of the 
Mobius algebra, and satisfies the recursion 

\B\ 

(21) = E 

i=l 

which is an immediate consequence of Eq. ([9]) and holds under the non-degeneracy conditions 
mentioned above. As introduced earlier, let A! be the partition that emerges from A by 
removing all singleton parts, which includes the limiting case 0 ' = 0 . 

Proposition 2. Let S be a finite set and V an indecomposable and irreducible lattice of 
partitions of S. For the generic choice of the recombination rates , the recursively defined 
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functions Mg(V(f7)) introduced above satisfy the relations 

for all A E Y(JJ) and all non-empty U C S. Here, A is any partition that emerges from A 
by the removal of an arbitrary number of singletons, including the case A = A!. 

Proof The claim is true for all A E P(5') with any S of cardinality |5| ^ 3 as a consequence 
of = p^{A,A) = 1 for all A E P(5) in this case, hence also for the sublattice V. The 

formula is also true for arbitrary [/ 7 ^ 0 when A contains no singletons or when A = A (both 
sides are equal) or when ^ = 0 together with A = A! = 0 (where both sides evaluate to 1 ). 

To continue the proof by induction in the number of elements of S respectively U, we may 
thus assume that A = Ty\J£ is the join of a partition H, which may or may not contain one or 
more singletons, with another (non-empty) partition £ that only contains singletons. With 
D := supp(T>) / 0 , we can then calculate 

\B\ |C| 

= E E qHc)11^^'{a\c,). 

i=l r)pV i=l 

Using the induction hypothesis and the structure of A, the last product can be simplified as 

|C| lU 

i=l j=l 

which allows to evaluate the inner sum over C as g^{H) — £<‘^(1) 6 {J-, Ij^))- This gives 

lU 

V^jF^l j=l 

= h^^hid) + q^{i\d) - /( I ) - 

where we used the recursion (j 2 in as well as the identity 'if^ifD) = 'if^{A) and the relation 
between if^d) and which both follow from the definition of the decay rates. Since 

7 ^ ip^iA) by assumption, our claim follows. □ 

Note that 1 ? coincides with the function ^ used after Corollary 0] to calculate the recombina¬ 
tion rate dependent terms in Table [TJ Inspecting this table again together with the recursion 
for 6 from Eq. m, one sees that a similar reduction as for i? must be valid for 9, too. Indeed, 
if A ^ B, one can replace the pair (A, B) by a reduced one that emerges via the removal 
of all common singletons in A and B. With O^{ 0 , 0 ) := 1, the reasoning of the proof of 
Proposition [21 applied to the recursion ([9]) for 9^, gives the following result. 
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Theorem 3. Let S and V be as in Proposition^ For the generic choice of the recombination 
rates , the functions 9^ £ A^{Y{U)) satisfy the relations 

for all A,B£ Y{U) with A ^ B and all 0 ^ U S. Here, B is any partition that emerges 
from B by the removal of an arbitrary number of singletons, including the case B = B'. □ 

This reduction property helps to better understand the values of the ^-coefficients and 
their recursive structure. Moreover, it also explains why and how the entries of the matrix in 
Eq. (11511 reappear within Tabled] for five sites, once for partitions that comprise the singleton 
{5} and once for those comprising {1}. 

Invoking Lemma [3| from the Appendix and observing that, for A ^ B, chn(A, B) and 
chn(p4.|g^pp^g,^,i3') have the same cardinality, it is clear that the reduction property for 6 
extends to its inverse p as follows. 

Corollary 6. Let the assumption be as in Proposition [2] and Theorem [S] and assume that 
the recombination rates are chosen such that 9^ G Aj^(V) is well-defined and invertible, with 
inverse . Then, one also has 

r^AB) = 

for all A,B G V with A ^ B. As before, B is any partition that emerges from B by the 
removal of an arbitrary number of singletons, including the case B = B'. □ 

Let us finally take a closer look at the structure of the probability vectors as a function 
of time, first in the simpler case that if = x-, with g^^^ := 6 {^)- 

Proposition 3. Let S be a finite set and consider the lattice I(S'). Moreover, assume that 
the recombination rates ^>(A) are chosen so that if = x on 1(5'). Then, the solution of 
Eq. ([6|) with initial condition aQ{A) = (5(A, 1) also solves the linear ODE system 

hW = - E Q{A,B)a^{B) with 

B>pA 

Q{A,B) = - MAC)x(C) = -gtotS{A,B) + ^ g{C) S{A,B AC), 

Ce[A,B] C^A 

where Q is {the transpose of) a Markov generator, with the diagonal entries Q{A, A) = —x{A) 
being the eigenvalues of Q. 

Proof. Under the assumptions, we know a^{A) = from Proposition [1] 

together with = Y1b)pCH{^) by Mobius inversion. Consequently, one has 

= -Y,L{A,C)x{C)e-^^^^^ = -Y,hiA,C)x{C)J2H{B) 

Of^A C>pA B-^C 

= E(- E h{A,C)x{C^H{B), 

B)pA^ CelA,B] ' 
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which proves the hrst claim and the first formula for Q. 

Clearly, Q{A, B) = 0 whenever A ^ B. If ^ ^ 0, one uses Eq. (fTOj) to get 

Q{AB) = - Y. ix{AC)Uot- Y 

C&[A,B] C&[A,B] ^ ^>6[C,1] 

=-e,,t5{AB)+ Ye{'D) Y 

'D)pA C€[A,V]n[A,B] 

Since [A, B] n [A, V] = [A, B AV], the last sum is bav] B AV) and the 

second formula for Q follows. 

All off-diagonal terms of Q are then clearly non-negative, while 

^ Q{A,B) = -Qtot+Y + E <^(A^AC) = 0 

.AGn(5) A^B C^A Cel(S) yl=^BAC 

because B AC is a unique partition in 1(5), so that the last sum is always 1. This shows the 
claimed Markov generator property. The diagonal entries are clear from the first formula for 
Q. They are the eigenvalues because Q is upper triangular. □ 

Note that the second expression for Q in Proposition [3] is slightly misleading in the sense 
that it suggests that all rates q{C) with B A C = A contribute to Q{A,B). However, the 
condition ip = X forces g to vanish for many partitions as |5| increases. In fact, as we shall 
see from Eq. (|2^ below, Q{A, H) > 0 is only possible if ^ ^ H and if A refines precisely one 
part of B. This is a somewhat surprising consequence of the underlying stochastic partitioning 
process discovered in O Sec. 6]. To explore this connection, one can now apply an analogous 
reasoning as in Proposition [3] to the general generic case, with the following result. 

Theorem 4. Let V he an indecomposable and irreducible lattice of partitions of a finite set 
and assume that the recombination rates q{A) are such that 9 € Ajj(V) is well-defined and 
invertible, with inverse rj. Then, the solution of the ODE system ([6]) with initial condition 
aQ(>I) = S[A, 1) also solves the linear ODE system 

aM) = with Q{A,B) = - Y 

B-^A Cg[AK] 

where Q is {the transpose of) a Markov generator. Moreover, Q is upper triangular with 
diagonal entries Q{A,A) = —ip{A), which are the eigenvalues of Q. 

Proof. The first claim follows from Eq. ([7|) in conjunction with Theorem [2j The calculation 
runs as in the proof of Proposition [3l now with 7y(C, B) a^{B) via the inversion 

formula for Eq. ©• Clearly, we have Q{A, B) = 0 whenever A ^ B, which implies that Q 
is upper triangular. The diagonal entries are then the eigenvalues, and their values are a 
consequence of r}{A,A) = 1 /6{A,A). 

Next, for fixed H € V, one finds the column sums 

YQiAB) = -Y E 0{A,C)iP{C)g{C,B) =Y^iC)g{C,B)Y(^i^^C)- 

4gV A^B C&[A,B] C^B 
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Since the last sum is <5(0, C) by [5l Prop. 8 (2)], which holds in this generality, one has 

= mriiQA) = 0 

from 1/5(0) = 0, as is necessary for (the transpose of) a Markov generator. Since Q is upper 
triangular, it remains to show that Q{A, B) ^ 0 whenever A B. 

To this end, we employ the underlying partitioning process from [U Sec. 6]. Our coefficients 
at{A) were identified as the first row of a Markov semigroup, namely as transition probabilities 
P^(l —>■ A). The corresponding Markov generator Q could be identified by also considering 
Pt{B A) for an arbitrary B €Y, leading to the so-called Kolmogorov backward equation. 
By an application of m Thm. 2.1.1], we can switch to the corresponding/orwarc? equation 
with the same generator, thus giving a linear ODE system for the quantities Pt{l —>• A) and 
hence a^{A) alone, which is nothing but the equation stated in the theorem, with Q = . 

Consequently, all off-diagonal elements of Q are non-negative. □ 

Since it is difficult to calculate the off-diagonal elements of Q from the formula in Theorem^ 
we employ [5l Thm. 3] to determine them. As mentioned earlier, Q{A,B) > 0 requires A ^ B 
with A refining precisely one part of B. So, we may assume this part to be B^, hence 
B = {B^, B 2 ,..., B^} with r = \B\ and A = {A ^,..., A^, B 2 ,..., B^}, with s ^ 2 and obvious 
meaning for r = 1, where |A| = s -|- r — 1 and A\q^ = {A ^,..., A^}. Then, we get 

(22) Q{A,B) = = E S{C), 

Cgv 

which clearly is always non-negative. One can now calculate that = V'('S)) 

line with Q having vanishing column sums. It is also instructive to calculate the matrix Q for 
the lattice V = II({1, 2, 3,4}) in both ways (via the formula from Theorem [4] and via Eq. (1221) 1 
to get an impression of the miraculous cancellations that happen along the way. The general 
verification of the equality of the two formulas is more difficult, and can partly be based on 
the reduction properties for 6 and rj, but also needs some lengthy calculations that we omit. 

Now, the Markov generator Q has a surprisingly simple form, and certainly no singularity 
of any kind. In fact, one can calculate it for V = P(<S'), where it holds for all values of 
the recombination rates g, and then restrict it to any sublattice, even to decomposable or 
reducible ones. The role of the somewhat tricky degeneracies appear in a new light then, too: 
Depending on the values of the parameters, the matrix Q will or will not be diagonalisable. In 
the former case, we are in the regime where the recursive solution with the ansatz from Eq. © 
works. In the latter case, we have to employ the well-known theory of linear ODE systems 
for Jordan blocks, compare [U Thm. 12.7] or [T71 Ch. 17.VII], thus giving the appearance of 
terms of the form T"e““* another and more transparent meaning. Eurther ramifications of 
this connection will be explored in a separate publication [3]. 
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Appendix: An inversion lemma for the incidence algebra 

Let L be a poset with ordering relation We assume L to be locally finite, which means 
that all intervals [A, B] in L are finite; compare [21 Chs. II and IV] for background. Now, let 
A^(L) be the corresponding incidence algebra (over any field K of characteristic 0). This 
is an associative A-algebra with the Kronecker function 5 as its (two-sided) identity. Any 
element of the algebra, a say, satisfies a{A, B) = 0 whenever A,B € L with A B, and the 
multiplication of two elements is defined by 

{a*fi){A,B) = a{A,C)f3{C,B). 

C&[A,B] 

An element a € A(L) is invertible if and only if a(A, A) A ^ for all A ^ L. In this case, its 
unique inverse a~^ (which is both the left and the right inverse, so a~^ * a = a* a~^ = d) is 
recursively given by a~^{A,A) = l/a{A,A) together with 

(23) a-\A,B) = ^^ Yl ^(AV)a-\V,B). 

It is sometimes useful to have a closed expression for a~^ that is not recursive in nature. 
One classic identity, see [161 Thm. 3.3.22], uses a summation over chains. If A 0, a chain 
from A to is any increasing sequence G C L of the form 

A = Qi < Q2 < • • • < Q^ = B, 

where ]G] = ^ is called the length of G (note that some authors use ]G] — 1 for the length). 
We use chn(A, B) to denote the set of all chains from A to B. When A ^ chn(A, B) = 0 
as usual. 


Lemma 3. If a ^ invertible, its inverse is given by 


a 


-\A,B) = Y 

Gechn(.4,B) 


-1 n!=i ^a{Gi,gi+i) 

rijJi ofiGjiGj) 


with the usual convention that an empty product is 1. 


Proof. The claim can be established by induction on the maximal chain length in chn(A, ;S). 
Clearly, the formula gives q;“^(A, A) = l/a(A, A), where chn(A, A) consists of G = {A} 
only. If A ^ B with the property that the maximal chain length between them is 2, which 
means that G = {A, B} is the only chain from A to B, the formula gives 


a-^{A,B) = 


a{A, B) 

a(A, A)a{B,B) 


which clearly agrees with Eq. (j23|) for this case. 
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Assume now that the formula holds for any pair with maximal chain length k ^ n, and 
consider A < B with maximal chain length n + 1. Then, Eq. (12311 gives the relation 


a-\A, B) 


E ct(^Aj C) 

c({A, A) 

A<C^B ^ ’ ' 


a-^{C,B). 


Now, by construction, the maximal chain length between any C and B in the sum is bounded 
by n, so that we can use the induction hypothesis as follows. 


a-\A, B) 


a(A,C) nl=l 

GecMC,B) ni='i 


= E (-1)'^' 

{^}UGGchn(^,B) 


«(Agi) nSi 

a{A,A) 


= E (-1)'^'-^ 

H£chn{A,B) 


nlS ^a{n,,n,+i) 


where the calculation really is a consequence of how each chain from A to B can uniquely be 
decomposed into one step from to C, with A < C ^ B, and a chain from C to □ 
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